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THE RANGE RATE ERROR DUE TO THE AVERAGING 
TECHNIQUES OF DOPPLER MEASUREMENTS 

by 

B. Kruger 
ABSTRACT 

The averaging error in slant range rate measurements, based on 
Doppler frequency integration over a finite time, has been analyzed. This 
error is of such a magnitude that a correction is required. For example, 
in a near circular earth orbit the error is 0.12 m/s for 370 km (200 nm) 
orbit altitude and0.5m/s for 185 km (100 nm) orbit altitude if an integra- 
tion time of 1 sec. is used. Means of correction, added noise due to cor- 
rection, and the residual error after correction are discussed. Special 
emphasis is given to the three point correction, which is shown to be 
adequate in most applications. In the above case the errors are reduced 
to 0.0001 m/s and 0.0024 m/s respectively. 
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THE RANGE RATE ERROR DUE TO THE AVERAGING 
TECHNIQUES OF DOPPLER MEASUREMENTS 

INTRODUCTION 

A finite time is always required for frequency measurements. The Doppler 
frequency shift can therefore not be measured instantaneously. Instead, the 
Doppler shift is integrated over a finite time T. This integral is proportional to 
a finite increment in range AR. The Doppler measurement, therefore, yields the 
average range rate R a during the time interval T 


R 

a 


T 



R dt 


AR 

T 


( 0 . 1 ) 


The actual range rate at time t* isR. 


The difference AR between R and R 

a 


AR = R - R 

a 


( 0 . 2 ) 


and suitable means for correction are discussed in this paper. The omission of 
the correction for AR will introduce an error in the interpretation of the measure- 
ments and this error is referred to as the averaging error in this paper. 


1. THE AVERAGING RANGE RATE ERROR 

The averaging slant range rate error AR can be calculated by expanding R 
in a Taylor series [1] around t = t Q 


* A 1 1 time in this paper refers to the spacecraft transponder rather than to the ground receiver in 
order to avoid corrections for propagation delay in this paper. The effects of propagation delay 
will be treated in a coming paper. 
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R(t 0 + At) .R 0 + R 0 At + Ir 0 At 2 +-|'R 0 At 3 + ... 


( 1 . 1 ) 




At = T /2 


At = T /2 


" T /2 


Figure 1— The Integration 
Interval 


*o + 1/2 


where 


R o - R' 


R 0 =» 


t = t 0 etc. 


If t 0 is the middle of integration interval T, then AT = T/2 as shown in Figure 1. 
The change in R during T is 


AR = R 



- R 



and applying equation (1.1) we obtain 


AR = R 0 + R 0 1 1 i F 0 




+ • • • 



+ • • 


or 


AR = R T + — R n T 3 
o 24 0 


and hence 


R 

a 


AR 

T 


R o + 24 R 0 T 2 + 


( 1 . 2 ) 


(1.3) 
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Using R a instead of R means that an error AR is introduced 


AR=R-R. = -UiS,T»- ... (1.4) 

It is seen from the above derivation, that all even derivatives R , d 4 R/dt 4 , 
etc. disappear in equation (1.4) if and only if, t Q is chosen in the middle of the 
integration interval T. Also the coefficient for the third derivative in equation 
(1.4) has a minimum for this choice of t Q . We therefore conclude that inter- 
preting R a as R in the middle of the integration interval minimizes AR. 

The averaging error is of such a magnitude that a correction is required 
in many cases. A complete correction is, of course, not possible. A practical 
limitation is that only one correction term is desirable, i.e., only the third 
derivative is used for correction and the fifth and higher order derivatives are 
neglected. A theoretical limitation is that all derivatives are contaminated with 
noise and a correction for ARwill therefore increase the noise in R. 

A practical important case is the circular orbit. A circular orbit will there- 
fore be used as a vehicle to demonstrate the magnitude of AR and of the residual 
error due to the neglect of the fifth and higher order derivatives. An overhead 
pass is chosen, because this is the worst case. 


2. THE AVERAGING ERROR IN A CIRCULAR ORBIT* 

A circular overhead orbit is shown in Figure 2. From Figure 2 we obtain 


R 2 = r 2 + (r + h) 2 -2r(r+h)cosi9 


( 2 . 1 ) 


(r + h) 2 = r 2 +R 2 + 2 r R sin 


( 2 . 2 ) 


s in 9 cose 


R r + h 


(2.3) 


From these equations d 3 R/di9 3 can be derived as shown in Appendix I. 


“This case is treated in Ref. (2), but due to analytical errors the results are not applicable. 


3 


Spacecraft 



i 

r = Radius of earth 
h = Altitude of spacecraft- 
R = Slant range 
6= Elevation angle 
6 = Central angle 


GODDARD SPACE FLIGHT CENTER 
SYSTEMS ANALYSIS OFFICE 
F«b 1965 


Figure 2-Orbiting Geometry 


In a circular orbit 6 is constant, therefore, 


Ad 3 

where R and d are derivatives with respect to time. AR is then, using the results 
from Appendix I, 


AR 




8 3 T 2 r 3 

8h! H) 



+ — 


r 


sin' 


cos e 


(2.4) 


4 




provided 


For e=0 


. 2 2 h 
s m e > > 


AR 3 


9 3 T 2 r 
24 


(2.5) 


With 


9- 


1/2 


(r + h ) 3 / 2 


( 2 . 6 ) 


where fi - the gravitational parameter, we obtain 


AR 


fJ /2 T 2 


8 h 2 r 3/2 (l + L 1 b 


2 r 


4 2 h . 2 , 
s i n e + — sit c] cos e 
r 


(2.7) 


provided sin 2 e ss 2h . 


I d 3 R/d0 3 | max is calculated in Appendix I, equation (A-ll). Using this value, 
we obtain 


ar 


2 M 3 / 2 T 2 

25 VS h 2 r 3/2 (l +3 — 


AR occurs for 

1 'max 


sin 6 = 



_h_ 
5 r 


( 2 . 8 ) 


(2.9) 
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Let us use the notation AR Q for the averaging error for a ’’standard" 
circular earth orbit with 


h Q = 1.85 km (= 100 nm) 
and a "standard" integration time 


The graph AR Q vs. the elevation angle e is shown in Figure 3. Inspecting Equa- 
tion (2.7), we find that 


AR^T 2 


and 


AR [except f o r cor rec t ion t e rms 

r 


For other h and T values than the "standard" values, we thus have 


AR 


f\j 

■ r X> 




( 2 . 10 ) 


Figure 3 shows | AR J = 0.5 m/s. I AR „| occurs at approximately 63° 

° 1 0 ‘max 1 0 1 max A 

elevation angle, which is in the region most important for tracking. 

A correction for AR will be required in many applications. If only the third 
derivative is used, as in Equation (1.4), a residual error will remain. The 
residual error is analyzed in the next chapter. 


3. THE RESIDUAL ERROR 

Terms containing d s R/dt s and higher derivatives have been neglected in 
derivation of equation (2.4) and (2.8). In order to estimate the residual error, 
we will find an upper bound for the magnitude of the term containing d s R/dt s . 
Including one more term in equation (1.4) yields 
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AR 0 (m/s) 



Computing AR for other h and T values, use GODDARD SPACE FLIGHT CENTER 

AR» (h 0 /h) 2 (T/T 0 ) Z AR 0 SYSTEMS ANALYSIS OFFICE 

Feb 1965 

Figure 3-Averaging error AR v s elevation angle £ for circular earth over 
head orbit with hg= 185 km (=100n. mi.) integration time Tq = 1 sec. 


J_ d^R 1 d s R 
2 ^ dt 3 1920 ^t 5 


(3.1) 


Using S 5 for the term containing d s R/dt s we obtain 


1 d s R 
5 = ' 1920 


(3.2) 
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In Appendix I, Equation (A- 11), it is shown that 


d s R 
de 5 


/ r 5 

< 23. 6 J_ 
h 4 


for an overhead pass. Hence 


Is, 


6 ' 


r 5 T 4 


81.4 h 4 r 5/2 [l + 


(3.3) 


or 


A 


5/2 j4 


81.4 h 4 r 5/2 (1 + 


(3.4) 


For an earth orbit with h = 185 km, the residual error is 


8 S I < 2.7-10- 4 T 4 m/s 


(3.5) 


The upper bounds |S 5 l max vs. T for overhead earth orbits are shown in Figure 4 
with h as parameter. It is seen that S 5 can be neglected for T = 1 sec. But for 
T= 3 sec., & 5 is more than 0.02 m/s for h = 185 km. Figure 4 thus clearly indi- 
cates the upper limit for the integration time. These values are, however, very 
optimistic because of the following reasons. 

The derivation of the residual error in this chapter was based on the assump- 
tion that d 3 R/dt 3 is known accurately. d 3 R/dt 3 is generally not known, but can 
be determined from a number of consecutive measurements. A minimum number 
of three measuring points is needed in order to determine d 3 R/dt 3 . Three con- 
secutive measuring points can therefore be used for the correction of the averag- 
ing error. d 3 R/dt 3 is, however, only determined to a certain accuracy by three 
(or any other finite number) measuring points and an increase in the residual 


8 


error can therefore be expected. In the next chapter the residual error for the 
three point correction is analyzed and found to be 9 times larger than the theo- 
retical residual error as given by equations (3.3) through (3.5) and as shown in 
Figure 4. The residual error for the three point correction is shown in Figure 7 
and 8. 



Integration time T (sec) GODDARD SPACE FLIGHT CENTER 

SYSTEMS ANALYSIS OFFICE 
Feb 1965 

Figure 4-Max. residual error |S- | „ after correction for d 3 R/dt 3 . I „ 

' d 'max 1 o 'm a x 

is due to neglecting d 5 R/dt 5 . d 3 R/dt 3 is assumed to be accurately known. 
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4. THREE POINT CORRECTION 


It was assumed in the previous chapters that the third derivative of R is 
known. This is generally not the case, but from three consecutive measurements 
R can be calculated and AR can be corrected for. The correction equation, the 
increase in noise, and the residual error for the three point correction are 
discussed in this chapter. 


The Correction Equation 


The equation for the three point correction can be obtained in the following 
way. It is assumed that three consecutive measurements are made with repeti- 
tion time T as shown in Figure 5. The measured average range rates are R al , 
R a2 , and R a3 respectively. The integration time is assumed to beT for all three 
measurements. This is not strictly true for the N -c ounting technique (3) in which 
T is a function of R a . During three consecutive measurements the variation in 
T is, however, small enough to be neglected. It is shown in Appendix n, Equation 
(A -33) that under the above conditions the range rate at time (t 0 + r) is 


R (t 0 + T) 


T 2 


24 T 


2 2T 


r t 


R 


l 2 T; 


T 2 T 2 


R 


1 2 T 2 T 2 


a 2 


T 2 


2 + ~ 2 t ~ 1 R a3 + S 4 + S s 


24 t 2 zij 2T 2 


(4.1) 


where S 4 and S $ are the residual errors due to the neglection of fourth and fifth 
order terms. It is also shown in Appendix II that 


5 


4 


I d ' R 

6 dt« 1 



— + r: X ) s. 

4Tj 2 T l) T, 


(4.2) 


10 


A 


and 


S, = 

5 24 


1 d5R - T 4 


dt 5 


T 2 7 T 4 
+ 

12 T 2 240 T 2 



( 4 . 3 ) 



Figure 5— Three point correction based on three consecutive measurements 


The Residual Error 

It is shown in Appendix I that for a circular overhead path 


and 


d 4 R 

d0 4 



d s R 

d(9 s 


^ r 5 

£ 23.6 I_ 
h 4 


For a circular orbit dd/dt is constant and therefore 


( 4 . 4 ) 


( 4 . 5 ) 


d n R _ d n R (d6 \ n 

dt n d 9 n ' dt ' * 4 * 6) 
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Using Equation (2.6) 


dd _ /J /2 

dt " (r + h) 3/2 


and Equations (4.2), (4.3), (4.4), (4.5), and (4.6), the upper bounds for the residual 
errors can be obtained. 



(4.7) 


M s/2 T 4 

| S | < 0.98 1 

j -^4 

A + t2N 

\ t2 

T 4 

h 4 r 5 / 2 (l + 7.5 £) 

12 T 2 240 T 4 

V 2T u 

/T 2 


which also can be written 





|S 4 I 





l‘.l 

; Is,! 

1 5 1 max 




For a circular earth orbit with p. = 3 

99 * 10 14 m 3 /s 2 , 

r = 6.38 • 

10 6 

m, Tj 


(4.8) 


and h = 185 km we obtain 


l*«l 


= 0.263 


, T 2 T 2 


r 

_ 

4T 2 T 2 

l l 


T i 


(4.9) 


K\ 


5 1 max 


0.0212 


12 T 2 240 T 4 



m/ s 


(4.10) 
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Equations (4.9) and (4.10) show the relationship between the residual errors 
and the times T, Tj, and r. The equations can also be shown graphically. Figure 
6 shows the residual maximum errors IsJ and IsJ for the case that R 
is referred to the time r. For convenient reference, the total measuring interval 
is shown at the bottom of Figure 6. It is seen, that if r = 0, i.e. R is referred to 
the middle of the measuring interval, then |5 4 | max = 0. For other r values 
I S 4 1 max increases and reaches a maximum at approximately r = ±0.6 T and 
decreases to zero at approximately r = ±1.1 T . The solid graphs are for the 
case T= Tj and the dotted for T = 1/2 Tj. The graphs are computed for a circular 
orbit and an overhead pass with an altitude of 185 km and T = 1 sec. 

It is seen from Figure 6 that R should be referred to the middle of the measur- 
ing interval, i.e. r = 0, because |sj max is zero. The two other zeros at r = 

±1.1 Tj are functions of T and the slopes of the curves are very steep in the 
neighborhood of these two zeros. Small variations in T or T will therefore 
produce relatively large residual errors. If, for instance, N-counting techniques 
[3] are employed, considerable variations in T have to be expected. 

We see also from Figure 6 that there is some residual error |S | for r = 0. 
Equation (4.8) shows that |§ s l ma![ is proportional to ( T,) 4 . The rapid increase 
of 1 8 | for increasing T is shown in Figure 7 for r = 0, T = T, and h as a 
parameter. For h = 185 km, |S s ! max is 4‘10 -2 m/s if the integration time 
T is 2 seconds. This residual error exceeds already the accuracy claimed for 
some range rate systems. Figure 8 shows |S 5 | max vs. h with Tj as parameter 
for the case T = T.. The rapid increase of |S_| for low values of h is evident 

1 1 j ma X 

from this graph. 


The R. M. S. Noise 


The total R. M. S. noise 77 after correction is obtained by adding the noise 
77 from each measurement in a R. M. S. sense. If we write Equation (4.1) in the 
form 


where 


R(At) = T 1 a, R „ + 8, + S s 



V 


(4.11) 


(4.12) 
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T = Tj (sec) 


GODDARD SPACE FLIGHT CENTER 
SYSTEMS ANALYSIS OFFICE 
Feb*"l965 


Figure 7-Three point correction of R(t n ) max. residual error §, 

U 1 5 'max 

5th and higher order terms. T is T, . 


due to neglecting 




Figure 8-Maximum residual error |Sj | mQx after three point correction vs. altitude h for the case 
of a circular orbit and an overhead pass. R is referred to the mid point of the measuring 
intervals (r= 0) and the integration time is equal to the repetition time of the measurements (T = T,). 
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we obtain the total noise 


V t = V '/z aj 


(4.13) 


assuming that each measurement has the same noise 77 . The ratio 77/77 is shown 
in Figure 6 for two values ofT: T = Tj and T = 1/2 Tj. It is seen that the noise 
increases 8 . 6 % for T = Tj and r = 0. The decrease of R. M. S. noise of 19% at 
r ~ 0.7 Tj can generally not be utilized because of the large residual error |S 4 |. 

5. CONCLUSIONS 

It has been shown that the averaging error for a circular earth orbit is 0.5 
m/s for 1 second integration time and 185 km orbit altitude. The averaging 
error is proportional to the square of the integration time and could therefore be 
reduced by shortening the integration time. The disadvantage of shortening the 
integration time is the increase in other errors such as quantization errors. 

The averaging error can also be reduced by applying corrections. A three point 
correction, based on three consecutive measurements, can reduce the residual 
error to below 0.01 m/s, provided that the corrected R is referred to a time cor- 
responding to the middle of the three measurements. The noise is shown not to 
increase more than 8 . 6 % in this case. The three point correction is therefore 
adequate in most applications. 

The time to which R is referred is rather critical in a three point correction. 
If the integration time varies or is not known, then a four or five point correction 
may have to be used. 
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APPENDIX I 


Derivation of R 111 

Differentiating equation (2.1) with respect to 9 yields 

RR 1 = r(r + h) sin 9 

(R 1 ) 2 + RR 11 = r(r + h) cos 6 

SR^ 11 + RR m = - r ( r + h) sin 6 

thus 



From equation (A-l) and (2.3) 


R 1 = r cos e 


is obtained. Solving for R n , using equation (2.1) and (A-2) 


r 2 sin 2 e + rh + — h 2 

R 11 = _ — 

R 2 


hence 


r 2 sin 2 e + rh + — h 2 
2 


R 2 




(A-l) 

(A-2) 


(A -3) 


(A -4) 


(A- 5) 


(A -6) 


(A -7) 
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Equation (2.2) can be solved for 1/R 


1 

R 


r sin e + ^r 2 sin 2 e + 2hr + h 2 
2 h r + h 2 


and 


1 _ 2 r 2 sin 2 e + 2 h r + h 2 + 2r sine A 5 " s in 2 e + 2 h r + h 2 
R 2 (2 hr + h 2 ) 2 


or 


• 2 h 
sm' e + _ 
1 r_ 

R 2 " 7 2 1 + h 

r 


(A-8) 


. , 2 zn 

for sin z e >> — 
r 

and thus 



4 2h 2 
sin e + — sin' e ' 


cos e 


(A- 9) 


valid for 


sin' 


e > > 


2h 

r 


For sin e = 0 
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we find 


*» 


R 


hi 


r 


The maximum of R 111 is found from Equation (A-9) 


R'"| = -iL iL U , 3h ' \ 

■" 25 /5 h ! V 2r/ 


for 


sin e 



h 

5r 


(A-10) 


(A-ll) 


(A-12) 


Estimation of |R IV | and |R V | 

Differentiating (A-3) two more times with respect to 6 yields 

3(R n ) 2 + 4R I R m + RR IV = - r(r + h) cos Q (A-13) 

10 R n R in + 5R r R IV + RR V = r(r + h) sin Q (A-14) 

R IV has extreme values for R v = 0. From Equation (A-l) and (A-ll) we obtain 
for R v = 0 

10 R 11 R m + 5R'R IV = RR 1 
Substituting R m from Equation (A-4) yields 


- 10 R 1 R 11 


+ 



+ 5 R 1 R 1 '' = RR 1 


(A- 15) 
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Equation (A -15) has two roots 


R 1 = 0 

or 

0=0 


and 


- 10 R 11 


+ 



+ 5 R IV 


R 


For R 1 = 0 we obtain from Equation (A- 13) 


piv _ -3(R n ) 2 - r 2 - rh 
R 


0=0 


and thus 


R 


IV I 


max 1 




+ r 


After elimination of R 11 between Equation (A- 13) and (A- 17) we obtain 


3 rIV _R 4 R i R 111 R 12 

2 10 “ R “ R 


or 


IR^I < JL 

8 
*f — 

R 1 R 111 

3 

+ — 

(R 1 ) 2 

1 15 

3 

R 

2 

R 


(A- 16) 


(A-17) 


(A-18) 


(A-19) 
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From 


R| > h, |R X | < r 


and 


|R ni | 5 0.86 J_ 
h 2 


we find that 


|R IV | <2.29 — 

' ma x 2 

h 3 


if lower order terms of h/r are neglected. |R IV | , is larger than |R IV | 

° 1 1 max 1 ° 1 1 max 2 

and again neglecting lower order terms we thus have 


R IV | 

1 max 


R IV | 

1 max 



(A-20) 


Because of the complex nature of the equations we will establish an upper bound 
for |R V | rather than deriving the exact expressions. From Equation (A- 1) and 
(A -14) we obtain 


R v I < IR 1 


R 11 R 111 

+ 5 

R 1 R IV 

R 


R 


Using Equation (A-5) and (A-ll) we find 


R v | < 23.6 ll 
h 4 


(A-21) 


neglecting the lower order terms of h/r. 


APPENDIX H 


Three Point Correction 


Expand R in a Taylor series around t 0 , including up to fifth order terms: 


E(it) =K 0 ^ 0 At + ^R 0 At 2 +7T R 0 At 3 + ^Rj'At 4 +^R 0 v At= + 


The derivatives are 


R(At) =R 0 + R 0 At + ^R 0 A t ^ + ^R-At 3 + — Rq At 4 + 


R(At) = R 0 +RoAt Rq V At 2 + iR v 0 At 3 


R(At) =R 0 + R”At + ± Rq At 2 


R IV (At) = R 3V + Rq At 


R V (At) = Rq 


For the first measuring interval, see Figure 4 we have 


r(- t 1+ I)-r(- T , -I) 


23 



and thus 


1 ••• T 2 1 ,, T 4 

R , = R(- T ) +— R (-T ) + — R v (- T ) — — 

ai v i y T 3 1 v 4 s t v 


5! 


16 


(A- 


and in the same manner 


R 


a 2 


R ,_Ir H 1 RV H 
0 ' 3! 0 4 5! 0 16 


(A- 
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(A- 


R(±Tj), R (iTi) and R v (±Tj) can be expressed in R q , R 0 etc. with the aid of 
Equation (A-22) through (A -24) by putting At = ±T r Solving for R Q , R Q and R ( 
we obtain 


K =R, —k T 2 R* - 1 


T 4 R" 


0 - a 2 4! -0 - 16 . 5 , “0 


(A- 


T 2 


y _ ^a3 “ ^al 1 

0 " 2 Tj ~JJ V 1 


T? + _ R IV 


(A- 


R 


. R , — 2R-+R i 

, a 3 a 2 ' a 1 1 


- 77 <- 2T i * t2 > rV o 


(A- 



Substitution into Equation (A- 22) yields 


R(At) = 


f At 2 

2T 2 

i 


At 

rr 


T 2 


24 T 2 
1 * 


• A At 2 2T 2 . 

|R . + 1 + i R 


a 1 


T 2 


24 T 2 

1 J 


a 2 


At 2 At 


k 2 T 


2 2 T, 


T 2 

24 T; 


R 


’a 3 


i 

6 





T 


3 

1 


1 

7 T 3 4 

T 2 

+ — 

(l + T 2 \ 

At 2 At 4 ') 

24 

(240 T 4 

12 T 2 

l 2T .y 

t . 2 * T .* J 


RV 0 T ! 4 


(A-33) 
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